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OPTIMAL CONFIGURATIONS OF LINES AND A STATISTICAL 

APPLICATION 

FRANgOIS BACHOC, MARTIN EHLER, AND MANUEL GRAF 


Abstract. Motivated by the construction of confidence intervals in statistics, 
we study optimal configurations of 2 d — 1 lines in real projective space MP® -1 . 
For small d, we determine line sets that numerically minimize a wide variety of 
potential functions among all configurations of 2 d — 1 lines through the origin. 
Numerical experiments verify that our findings enable to assess efficiently the 
tightness of a bound arising from the statistical literature. 


1. Introduction 

Motivated by a question arising from statistics related to the construction of 
uniformly valid post-selection confidence intervals ®E], whose details we present 
below, we aim at computing 2 d ~ 1 evenly-spaced lines through the origin in R d . The 
terminology “evenly-spaced” is loose and, indeed, there are several mathematical 
formulations that make sense but often lead to different types of line sets. Usually, 
one considers a potential energy with a “repulsive” pairwise interaction kernel, and 
a line set may be called evenly-spaced if it is optimal with respect to this energy. 
However, note that optimal line sets may differ among different potential energies. 
For very particular choices of the number of lines with respect to the ambient di¬ 
mension d, see 0, optimal lines coincide for a large class of monotonic pairwise 
interaction kernels considered in [7] and are therefore called universally optimal. 
However, universal optimality is a rare event in the sense that only very few con¬ 
figurations can exist and even fewer are actually proven to be universally optimal. 
Indeed, in general it is extremely difficult to prove that a certain configuration is 
universally optimal, cf. dill 13- 

In the present paper, we consider configurations of 2 d — 1 lines in R d for d = 
2,..., 6 that minimize (not necessarily simultaneously) three types of potential en¬ 
ergies associated to the distance-, the Riesz-1-, and the log-kernels. Additionally, 
we compare these minimizers to the corresponding best packings of lines found in 
m, which can be seen as limiting cases of potential energy minimizers. In dimen¬ 
sion d = 2 , 3 , all the three minimizing configurations of 3 and 7 lines, respectively, 
are known to coincide with the corresponding best packings of lines, by virtue of 
their universal optimality [§]. Unfortunately, for d = 4 , 5 our numerical experiments 
suggest that there are no universally optimal line sets. Surprisingly, for d = 6 there 
seems to exist a universally optimal configuration of 63 lines, which we identified as 
the best packing configuration provided in m • However, this particular configura¬ 
tion, which can be composed by the 36 lines going through the vertices and the 27 
lines going through the centers of the 5-faces of the I 22 polytope (also called the Eq 
polytope), cannot be proven to be universally optimal with the so-call sharpness 
condition introduced in [7], so that we state its universal optimality as a conjecture. 
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Let us now present the statistical application. In the context of the construc¬ 
tion of valid post-model-selection confidence intervals, 00 have proposed new 
confidence intervals that are derived from evaluating a special statistical potential 
function of at most 2 d — 1 many lines in R d . It is desirable to find the maximum 
of this function, when the lines are subject to certain restrictions Eiia. Due to 
the inherent complexity of these restrictions and of the statistical potential itself, 
direct optimization approaches seem hopeless. However, a standard upper bound 
is available for this maximum derived by two consecutive inequalities, cf. <®> 
and 0 in Section [2j It is known that the second inequality is an equality for 
d = 2 j| Lemma A.4] and that it is tight as d —> oo 0. In the present paper, 
we shall complete this picture for other small values of d by simply evaluating the 
statistical potential at sets of 2 d — 1 evenly-spaced lines. These evaluations are very 
close to the upper bound, which demonstrate the tightness of the second inequality. 
Moreover, for universally optimal line sets, the gap is significantly smaller, which 
underlines that special property. 

The outline is as follows. In Section [2j we fix notation and present the statistical 
potential that motivates us to search for 2 d — 1 evenly-spaced lines in R d . The 
potential energies in projective space are introduced in Section [3j where we also 
provide the definitions of the distance-, the Riesz-s-, and the log-energy m as well 
as the notion of universal optimal configurations of lines 0. For each dimension 
d = 2,... ,6 we provide in Section [4] one of the numerically found minimizers of the 
distance-, the Riesz-1-, or the log-energy. Section 0 is dedicated to the statistical 
application, in which we compare the performance of the evenly-spaced lines with 
a naive Monte Carlo optimization of the statistical potential. 


2. Notation and motivation 

Let S d_1 denote the unit sphere and RP d_1 the projective space (the set of lines 
through the origin) of R d . Any u £ S d_1 defines a line i £ RP d_1 by £ = uR, and its 
antipodal counterpart —u yields the same line. Throughout the entire manuscript, 
we shall fix TV := N(d) := 2 d — 1. The set of all sets of at most TV lines is denoted 
by 2 ?<jv := T^<N{d) := {L C RP d_1 , < TV} and the set of all sets of exactly TV 
lines is denoted by Vjy := {L C RP d_1 , = TV}. 

In order to construct valid confidence intervals in statistical model selection, the 
authors in El El consider a function 


fd,r,ct '■ T^<N — t R+, 


where a £ (0,1) and r £ N* are fixed parameters. The value of fd,r,a(L), for 
L £ D< jv, is defined as the unique K > 0 such that 


(1) 


E, 




K 2 


'd- max (u.V) 2 ' 

■uR=(eL 


= 1 — a, 


where is the cumulative distribution function of the F-distribution with pa¬ 
rameters d and r, and V is a uniformly distributed random vector on § d_1 . In the 
statistical context, for fixed values of a and r, it is desirable to find the maximum 
of fd,r,a on a certain subset V C P< n, cf. 00, i.e., one aims to determine 

(2) SUp fd,r,a (-^) > 

Lev 
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where each set of lines in V is derived from some statistical data set but where 
D depends only on d, cf. [3 Equations (5.2) and (5.3) and Section 4.10] and [3] 
Equation (6)]. Note that [3] and [5] yield two different subsets V, but this difference 
is of no consequence in the sequel. Considering the supremum ([2]) is beneficial 
because this supremum is data-independent and can be tabulated, for fixed r and 
a, once and for all. 

Exactly determining © is difficult if not computationally infeasible since the set 
T> is defined in [3t or [3 in an intricate manner, that has so far made it impossible 
to determine Q theoretically. Moreover, the values f d , r ,a{L) can usually only be 
approximated by Monte Carlo methods, as in Algorithm |TJ which we write only for 
the case of sets of N lines, for concision. Note, that Algorithm fl is also used in [3]. 

In [3J [5] , an upper bound K (d, r, a) has been proposed for ( 21) that can easily 
be determined numerically, and we refer to Proposition 2.3 and Algorithm 4.3 in 
(3| for its definition and computation. This bound satisfies the following sequential 
inequalities: 

(3) sup fd,r,a(L) < sup fd,r,a(L) 

LeD L£V n 

(4) <K(d,r,a), 


where (|3| is due to T> C 2?<iv and to the fact that K in |l]) is increased if additional 
lines are added, so that sup ieX , SJV f d , r ,a{L) = sup LgX)jv f d , r ,a(L). The inequality 
Q is derived from a union bound, cf. Equation (A.17) in [5|. 

As noted in [5J and [3] Remark 2.10], the inequality Q is tight for large values 


of d , i.e., 

(5) 


SU PLgD w fd,r , a ( L ) 
K(d , r, a) 


when d tends to infinity, cf. [5] proof of Theorem 6.3]. For d = 2, we even have 
su Plg:d n fd,r,a{L) = K(d,r,a), see [3J Lemma A.4], 

In the present paper, we shall assess whether the upper bound K(d,r,a ) in Q 
for sup LGT , N fd, r ,a{L ) is also tight for small values of d. 

To evaluate the quality of the inequality Q, we must approximate 

(6) sup f d ,r,a(L) 

LeT>n 


to determine its difference to K(d,r,a). However, exactly computing |6]) is also 
difficult, if not numerically infeasible, and we shall rather aim to derive good lower 
bounds. In fact, we shall verify that evaluating f d , r ,a on one or f gw candidates 
of evenly-spaced lines L = {£\,... ,1^} C KP^ 1 yields better lower bounds on 
Q than several Monte Carlo attempts. Our numerical results on this issue are 
presented in Section [5] 

Remark 2.1. The quantity K 2 in [3] corresponds to ([2| where the supremum holds 
over a certain subset T> that depends on other quantities beside d (it is data de¬ 
pendent). In the case of K 2 , there exists another upper-bound K 3 in [3] , which 
is smaller than K(d,r,a), incomparable with sup igI3 fd,r,a(L ) and convenient to 
compute. Nevertheless, the context of this paper, where the subset T> is data inde¬ 
pendent, is unrelated to K 3 , so that K(d,r,a) is the only available upper-bound for 
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3. Potential functions in real projective space 

We shall now consider families of potential energies whose minimization can 
provide us with rather evenly-spaced lines. We recall that the chordal distance d c 
between two lines £i, £j £ RP d_1 is given by 

(7) dUliJj) = i ||muj - UjuJ || 2 , 

where £i = u,R and £j = Uj R with Ui, Uj £ S d_1 . Let / : (0,1] —> R be a decreasing 
continuous function. For fixed N and d, we aim to minimize the potential energy 

T-’f : 2?jv —t R 

i.e., we aim to find L £ T>n such that 

V f (L)<Vf(L), for all L £ V N . 

We shall explicitly consider three types of pairwise interaction kernels, 

(A) the distance-energy, 

fi(t) = -Vt, 

(B) the Riesz-s-energy , 

hit) = t~ s/2 , for s > 0, 

(C) the log-energy, 

h{t) = - log (t), 

and we refer to [si si nu na nu cm na and references therein, for investigations on 
their minimizers and asymptotic results when N tends to infinity. 

We call a function / : (0,1] —>• R completely monotonic if it is infinitely often 

differentiable and (—1 ) fe /( fe ) > 0, for k = 1,2,_ Note that the above fi, hi h 

satisfy this property. The following definition is borrowed from [7|. 

Definition 3.1. We call L £ T>n universally optimal if it minimizes Vf among 
T>n, for all completely monotonic functions /. 

Additionally, it seems natural to consider configurations solving the packing 
problem of N lines in RP d_1 , i.e., sets of lines which maximize the minimal chordal 
distance, i.e., find {i\,... , £n} £ R’n such that 

min d c (£i, £j) > min d c (£i,£j), for all {£ 1; ... ,£n} £ T^n, 

cf. [10] • Note that the packing problem corresponds to the limit case of the Riesz- 
s-potential when s tends to infinity. In particular, if the set {£\,... ,£n} £ Thv is 
universally optimal, it also solves the packing problem. 

For packings in Grassmannians beyond the projective space, we refer to PM, 
see also [2]. 
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4. 2 d — 1 EVENLY-SPACED LINES IN SMALL DIMENSIONS 

Recall that we are interested in minimizers of the distance-, Riesz-1-, and the log- 
energy for N = 2 d — 1 with d = 2,.. ., 6. Note, that for each potential energy there 
exist lots of local minimizers and, unless all such minimizers have been determined, 
there is no general statement that one has already found a global minimize!'. 

Here, we repeatedly apply a local optimization procedure initialized by randomly 
chosen starting points (in our case the nonlinear CG method described in mm), 
and we are convinced that the line configurations we present are actually minimizers 
of the corresponding energies. 

In general the minimizers of each of the three potential energies are different. 
However, since the different configurations perform rather equally well for the sta¬ 
tistical potential we will provide for each dimension only one selected minimizer 
explicitly. 

4.1. Universally optimal lines in R 2 . It is known that three equiangular lines 
in R 2 are universally optimal, cf. [7J, hence all minimizers coincide with that best 
packing configuration. For instance, such lines {^fc}| = i are given by the vectors 
Uk = (cos(|/c7r),sin(|fc7r)) T , with Ik — UfcR, for k = 1,2,3. Note that this config¬ 
uration is highly symmetric and generated by a single orbit of its symmetry group 
Z?3, which is the dihedral group of order |£>31 =6. 

Remark 4.1. It is known that those {It,^ 2 ^ 3 } maximize f2,r,a and, as mentioned 
before, the maximum indeed coincides with K(d,r, a), cf. [2i Lemma A.4]. 

4.2. Universally optimal lines in R 3 . It has been proven recently that there 
exist 7 lines that are universally optimal in RP 2 , see j9j, hence all minimizers 
coincide with that best packing configuration. In such a configuration, 4 lines are 
going through the vertices and 3 through the centers of the faces of a cube centered 
at the origin, see Figure [l] If the cube’s edges are aligned with the coordinate axis, 
then the lines {£k = WfcR}£ = i are given by the vectors 

Ul = (l,l,l) T , it 2 = (1, —1,1) T , U 3 = (—1,1,1) T , 114 = (—1, —1,1) T 

u 5 = (A 0,0) T , uq = (0,1,0) T , u-j = (0,0,1) T . 

Note that this configuration is also highly symmetric and is generated by only 
2 orbits of its symmetry group Oh, which is the full octahedral group of order 
\Oh\ = 48. 


4.3. Optimal lines in R 4 . Our numerical computations provide us with strong 
evidence that there are no universally optimal configurations of 15 lines in RP 4 . 
We found a very symmetric configuration L = d which seems to minimize 

the log-energy and the Riesz-1-energy simultaneously. It is more symmetric than 
our numerical minimizer of the distance-energy in the sense that it is composed by 
fewer group orbits. Moreover one can simply check that it is a stationary point of 
any Riesz-s-energy, s > 0. However, it does not solve the best packing problem 
(the limiting case s —> 00 ), as it has a slightly smaller minimal distance than the 
configuration found in [TOlj . 

The symmetry group Q of L as a subgroup of the orthogonal matrices 0(4) C 
R 4x4 acts naturally by left multiplication on RP 3 . The group has order \Q\ = 144 
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Z, 



Figure 1. 7 universally optimal lines in R 3 . 


and is generated by the following matrices 


l l 

0 

0 



(° 

0 

0 




0 

0 


0 

1 

0 

0 

, G 2 = 

0 

0 

1 

0 

, G 3 = 

0 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

1 

0 

V° 

0 

0 

V 


V 

0 

0 

(V 



0 

0 

V 



f - 1 

- V 3 

0 



(~ l r- 

-V3 

0 

0 \ 

II 

to | 

V 3 

0 

-1 

0 

0 

1 

0 

0 

, c 5 = l 

V3 

0 

-1 

0 

0 

-1 

0 

-73 


\0 

0 

0 

V 


l 0 

0 

73 

-1/ 


The 15 lines L = Li U L 2 are composed by the two orbits 


L 1 = {uR\u = G- (1,0,0,0) T , G e g}, 
L 2 = {uR\u = G ■ (0,1,1,0) T , G g Q} 


of cardinality 6 and 9, respectively. See Figure [2] for a visualization. 


4.4. Optimal lines in R 5 . We numerically derived different minimizers for each 
of the three potential energies. They all have a nontrivial symmetry group, and we 
want to provide the most symmetric configuration of lines, which is found for the 
Riesz-l-energy. In that case the 31 lines L = {Zk}\ Li are composed by 8 orbits of 
its symmetry group. More precisely, the symmetry group Q has order \Q\ = 12 and 
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Figure 2. 15 lines in R 4 visualized in the unit ball in R 3 . Blue 
points correspond to orbit L±, yellow correspond to L 2 . For visu¬ 
alization, we consider the intersection of lines in R 4 with the upper 
hemisphere of the unit sphere. Hence, lines become points. Now, 
we apply the stereographic projection to map the upper hemisphere 
in R 4 into the full ball in R 3 , which we can plot. In general, a line 
in R 4 reduces to a single point in the ball but lines that intersect 
the equator have two intersection points and we plot both. 


is generated by the following three matrices 


0 

1 0 

0 0 



0 

-1 

0 

0 

0 

0 

0 1 

0 0 

5 

g 2 = 

0 

0 

1 

0 

0 

0 

0 0 

1 0 



0 

0 

0 

-1 

0 

\o 

0 0 

0 1 



\0 

0 

0 

0 

V 

( 2 

0 

0 

0 

0 \ 






0 

-1 

-Vs 

0 

0 






0 

\/3 

-1 

0 

0 






0 

0 

0 

-1 

-V3 






\0 

0 

0 

v/3 

-1/ 
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The 31 lines L = 
Li = {zzR | 
L 2 = {zzR | 

L 3 = {zzR | 
L 4 = {zzR | 
L 5 = {zzR | 
Lq = {zzR | 
L 7 = {zzR | 
L 8 = {zzR | 


Ui = i Li are then given by the eight orbits 
zz = G-(1,0,0,0,0) t ,Gg£}, 
zz = G • (0,1,0,0,0) T , G & G}, 

u = G • (0, a\, 0, — \J\ — af, 0) T , G G G}, 

u = G ■ (&!, 0 , 62 , 0 , ^l~bl-bl) T ,G € G}, 

u = G■ (ci,0,-c 2 ,0,-yi^cf^c|) T ,G G Q}, 

u = G ■ {d u d 2 , -d 3 , di, ^l~d\-dl-dl~dl) T , G G G}, 

u = G ■ (ei, e 2 , e 3 , -e 4 , 1 - e? - e| - e| - e|) T , G G 1?}, 

u = G • (fi, -f 2 J 3 , U,y/l- fl - f! - fI - f!) r ,G G G}, 


of size 1, 3, 3, 3, 3, 6 , 6 and 6 , respectively, where the constants can be computed 
to arbitrary precision by numerical minimization 


ai = 0.1386569..., 
61 = 0.6107676..., 
ci = 0.6319241..., 
d x = 0.0959289..., 
d = 0.3556342..., 
fi = 0.5842996..., 


b 2 = 0.2652528..., 
c 2 = 0.6489064..., 
d 2 = 0.6048195..., 
e 2 = 0.0311526..., 
f 2 = 0.4841533..., 


d 3 = 0.3121361..., 
e 3 = 0.7210732..., 
f 3 = 0.3040229..., 


d 4 = 0.6715440..., 
e 4 = 0.4584122..., 
f A = 0.5287425... 


4.5. Optimal lines in R 6 . Our numerical investigations provide us with evidence 
that there exists an arrangement of 63 lines which is universally optimal in RP 5 and 
thus coincides with the best packing solution presented in m- This configuration 
L = {ik}t=i is given by the 36 lines going through the vertices and the 27 lines 
going through the centers of the 5-faces of the I 22 polytope, cf. Ill], also known as 
the E 6 polytope. The symmetry group G of that polytope has order \G\ = 103680, 
which contains the Coxeter group E 6 as an index 2 subgroup. The high order of 
the symmetry group G shows the remarkably high symmetry of this configuration 
of lines. Note that the group G can be generated by 2 matrices 

/I 0 0 0 


Gi = 


( 


\ 


0 0 


V8 

1 

V2 

_3^ 

y/8 


V 8 

VI 


V 2 

0 


v's 

v'f 


0 \ 

vT' 

0 

v/I 

0 


G 2 = 


0 

0 

0 

0 

Vo 


0 

1 

0 

0 

0 


0 

0 

0 

0 


0 
0 
0 
0 

1 73 

2 2 

\/3 1 

2 2 / 


so that the 63 lines L = L x U L 2 are explicitly constructed by the two orbits 
L x = {zzR | u = G ■ (1,0,0,0,0,0) T , G G E 6 }, 

L 2 = {uR\u = G■ ( 0 , 0 , 0 , 0 , Vs, 1 ) T , G G E 6 } 

of cardinality 36 and 27, respectively. Alternatively, the 63 lines can be derived 
from the minimal vectors of the E 6 lattice and its dual lattice E£, see, for instance, 
G. Nebe’s and N. Sloane’s website [ 1 |. 
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Note that L\ itself satisfies the sharpness condition in the sense of [7] and hence 
is universally optimal in MP 5 . The 63 lines are not sharp, but we conjecture that 
they are universally optimal in MP 5 . Moreover, the corresponding configurations 
of 126 points on the sphere is stationary for any completely monotonic potential 
function of the squared distance of 126 points on § 5 , cf. ( 2 j. 


5. Applications to the statistical potential 

5.1. Two approaches. Now that we have some configurations of N = 2 d — 1 
evenly-spaced lines in projective space in hand, we can apply them to derive lower 
bounds on 

(9) sup fd,r,a{{h,---jN}) 

{^i ,...,£n}$=T>n 

as was our aim stated in Section [2j We shall compare evaluating fd.r.a at our 
evenly-spaced lines with a standard Monte Carlo optimization method. Note that 
fd, r ,a cannot be evaluated directly, and we apply Algorithm [T] to do so, which is 
the same algorithm as [3, Algorithm 4.1], and which is also used by the authors of 
[5]. To summarize, we shall compare the following two methods: 

(i) evenly-spaced lines: 

This method simply consists in applying Algorithm[l]with {$i,£^} being 
one of the optimal line sets derived in Section [4] 

(ii) naive Monte Carlo: 

One aims to maximize fd, r ,a by Monte Carlo optimization, given in Algo- 
rithm[4j which first calls the random line generator in Algorithm [2] and then 
applies Algorithm [3j which itself repeatedly calls Algorithm [T| This type 
of Monte Carlo optimization has already been used in [3] Section 5] (in a 
context unrelated to this paper and the evenly-spaced line method (i)). 

Our intention is to check that (i) indeed outperforms (ii). From a computational 
point of view, (i) has a clear advantage: the set of evenly-spaced lines is obtained 
numerically once and for all, and can be used for any value of r and a. In contrast, 
one needs to repeat the naive Monte Carlo optimization for each values of r and a 
under consideration. 


Algorithm 1 Evaluation of fd,r,a at {ti,. .. ,£n} 

Input: lines {t\, ...,£n} £ Thv and some parameter I £ N 
Output: K approximating fd, r ,a({£\, £n})- 

1: generate independent uniformly distributed random vectors {V );}( =1 C 
2 : for each i = 1,..., I, calculate Ci = max.,- = i j ... i jv{'Uj, Vi) 2 , where tj = itjM. 
3: determine K that solves 

( 10 ) - Fd,r=l — a. 

2—1 1 

4: return K 
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Algorithm 2 random line generator 

Input: Ni and a set of probability distributions {Vj }jL-\ on S d_1 . 
Output: D\ C T>n with = N\. 

1 : initialize D\ = 0 
2: for i = 1,..., N\ do 

3 : for j = 1.... . N, generate a Uj -distributed random vector Uj £ S d_1 

4: insert the corresponding line set {iqR,..., uatM} into D\. 

5: end for 
6: return D\ 


Algorithm 3 find max 

Input: D\, parameters N 2 ,Ii,I 2 ,I 3 £ N with N 2 < #D\. 

Output: approximation of max^q fd,r,a{£ l, • • • ,£n)- 
1 : evaluate fd, r ,a at all {£ lt ..., £n} £ with I = I\ in Algorithm |Tj 
2 : keep the N 2 sets D 2 C D\ that yield the largest evaluations in step 1. 

3: reevaluate fd, r ,a at all {ii, ...,£n} £ D 2 with I = I 2 in Algorithm]!] 

4: keep the sets (f*, ...,£* N } £ D 2 that yields the largest evaluation in step 3. 

5: return K derived from the reevaluation of fd,r,a{{£*, Vl) Algorithm]!] 
with I = I 3 . 


Algorithm 4 naive Monte Carlo 

Input: parameters Ni, N 2 , I\, I 2 ,13 £ N with N 2 < N±. 

Output: naive approximation of supp^ fd,r,a({£ 1 , • • • 

1 : generate Di C T>n by calling Algorithm [2] with Ni and the uniform distribution 
on S d_1 . 

2 : return K derived from calling Algorithm [3] with D\ : N 2 , I\, I 2l I 3 . 


5.2. Numerical experiments. We compare the two methods presented in Section 
5.1[ for d, = 3,4, 5, 6. Two values of r = 20,60 are chosen as representative of 


practical statistical situations, and a = 0.05 is chosen since this yields the classical 
95% confidence. 

In order to further investigate on the approach (i), we shall also consider local 
searches in proximity to the evenly-spaced lines: 

(a) local evenly-spaced lines: 


We search for the maximum of fd, r ,a locally around {fi,..., given as in 
(i). Indeed, we call Algorithm [ 2 ] using N projected Gaussian distributions 
with isotropic variance a 2 = 0 . 1 2 and mean vectors {u*}f =1 , respectively. 
With the resulting set V 1 , we apply Algorithm [3] 


(b) very local evenly-spaced lines: 

As in (a) but we search even more locally by choosing <r 2 = 0.01 2 . 

In the numerical experiments, the parameters for (i), (ii), (a), and (b) are chosen 
by / = 20 000 000, W = 200 000, N 2 = 20 000, h = 10 000, I 2 = 200 000, and 
I 3 = 20 000 000. In Figure [3] we report, for each of the configurations of d, r (8 
configurations in total), the ratios K/K{d, r, a), where K takes four values obtained 
by the methods (i), (ii), (a), and (b), and where K(d,r,a ) is the upper bound in 
<§• For the methods (i), (a), and (b), for d = 3,6, a unique set of 2 d — 1 lines 
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is under consideration, which minimizes all the potential functions of Section [3] 
For d = 4, 5 different sets minimize different potentials, but give approximately the 
same values for fd,r,a (the lines corresponding to the packing problem slightly lag 
behind the others). Hence, in Figure [3] below, we only report the results for the 
minimizer of the Riesz-1-potential, for concision. The ratios K/K(d,r,a) enable 
us to compare the two methods, (i) and (ii), while those obtained from (a) and (b) 
address the local optimality of approach (i). 

In Figure [3j for the configurations of d and r under consideration, the method 

(i) using evenly-spaced lines provides a better maximization of fd, r ,a over than 

(ii) . Hence, it is beneficial both from a computational and performance point of 
view. 

We also observe in Figure [3] that the values of I\ obtained from (a) are below 
those obtained from (b). Furthermore, although not perceptible in the figure, we 
always have that either K is smaller in (b) than in (i), or the two values cannot 
be distinguished, because of the very small but positive variance of Algorithm [Tj 
This is numerical indication that the sets of evenly-spaced lines are local maximiz¬ 
ers for fd,r,a over T>n- At least in the cases d = 3,6, since we are then dealing 
with universally optimal lines, it seems plausible that we even obtained the global 
maximizers. 

Note that, in light of Figure [3| K(d,r,a) is a tight upper bound in Q for 
d = 3, ..., 6 , with a difference of less than 0.5%. This result is of complementary 
nature with the tightness for large d , see ([5]), derived in j5. : , proof of Theorem 6.3]. 

We can conclude that it would not be very beneficial to aim at improving the 
union bound Q. Instead, one may want to study the inequality § more closely. 
Given the complexity of the sets D considered in [3] or [5] , however, this may require 
other tools and may turn out to be an extremely challenging task going beyond the 
scope of the present paper. 

We have a smaller window of possible values for in dimensions in which 
universally optimal lines exist, or are conjectured to exist, with the predefined 
cardinality N = 2 d — 1, see Figure [ 3 ] where d = 3,6 yield higher ratios for (i) than 
d = 4, 5. In this sense, our statistical application indicates that the concept or 
property of universal optimality is indeed beneficial. 
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